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and, therefore, 

11 ifh. fy)= n (<P'y. P-'y') — n (<P'y'> J / — f'v) 

= n{<p' y ,,n>y) + Il(<p'y,-iy + fi'y), (16) 

or, writing for brevity fi y for the co-amplitude of ft', 

17 ifh, fy ) = n if'y'.P-'y') — n (<p'y, p>y) ; fly + p\> = ->y. (16') 

This I call Somoff' s theorem. Somoff, however, makes no use of it in 
this form. He employs a relation which is obtained from (16) by expressing 
II{(p'y' t Jyi — fi.'y>) in terms of Il(<p'y',fi'y') by (6). We have thus 

n {f y,,y ) = 2 n {f 'yy y) - fVrs^^f + ¥ ^yT^rPr 

nit , \ ■ , 1 z 1 + 7" sin usin <p . . 

= 2lJ (<p'y,/x'y')—r<fy sm ft + \l n ' . r . r . (17) 

KT ' ' T ' r i 1 — y sin // sin <p v ' ' 

[to be continued.] 



A SIMPLE DISCUSSION OF LOGARITHMIC ERRORS. 
By Prof. H. A. Howe, Denver, Col. 

[CONTINUED FROM VOL. II, PAGE 43.] 
PROBLEM XI. 

To find the probable error of an interpolated logarithmic trigonometric func- 
tion. 

Tables of logarithmic trigonometric functions may be divided into three 
classes. The first class embraces those in which the tabular difference is to be 
multiplied by a decimal of one, two, three, or, occasionally, four places. In this 
class we find 5, 6, and 7-place tables in which the functions are given for each 
10 seconds or for each second, those tables in which the functions are given for 
each minute, and are interpolated for decimal parts of a minute, the 4-place tables 
in which the functions are given for each 10 minutes, and tables in which the dif- 
ferences of the successive values of the argument are tenths or hundredths of a 
degree. 

The second class includes those tables in which the functions are given for 
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each minute, and the proportional parts for 6, 7, 8, 9, io, 20, 30, 40, and 50 sec- 
onds. Gauss's or Hoiiel's 5-place tables are such. 

In the third class are those 5 and 6-place tables in which the functions are 
given for each minute, and the variations of the functions for each second, these va- 
riations being carried to two decimal places. In Olney's or Davies's 6 place tables 
these variations, though printed directly opposite the functions, should be placed 
half a line lower, since they have been so computed that they show the rate 
at which the function is varying at each half minute. In Schlomilch's 5-place 
tables the variations have been found by dividing the differences between suc- 
cessive tabular functions by 60. 

Class I. The discussion of this class has already virtually been given in 
Problem I. If x be taken as 0.01, 0.02, . . . 0.99 the average probable value 
of xc' + (1 —x)c becomes 0.810 X 0.25, and no appreciable change will be 
made if we have x equal to o.OOi, 0.002, . . . 0.999. But taking into account 
the fact that x may be 0.0, in which case the probable error is 0.25, we find that, 
when x varies from 0.1 to 0.9, the average probable error is 0.3 1. When x va- 
ries from 0.01 to 0.99 and from 0.001 to 0.999, the average probable error is 0.32. 

Class II. In Hoiiel's 5-place table, the functions in the tables of proportional 
parts are given to the nearest integer only. If the numbe* of seconds has two 
figures and does not end in a cipher, two interpolations from the table of propor- 
tional parts must be taken, and the probable error of the sum of these is j/2 .(0.25), 
or 0.35. If n denote the number of seconds and x= n/60, the error of the in- 
terpolated function would be (aside from the source of error just mentioned), 
xc' -\- (1 — x)c. Since the probable values of ^'and c are each 0.25, we have the 
following formula for the probable error in this particular case : — 

Error = ^/\x % + (1 — xf + 2] . (0.25). 

When but a single number is taken from the table of proportional parts the 
proper formula is 

Error = v 7 ^ 2 + (1 — xf + 1] . (0.25). 

Substituting the values of x, which range from 1/60 to 59/60 in the preceding 
formulae, we find that the average of the probable errors is 0.39, the minimum 
being 0.31 and the maximum 0.43. 

But, if the functions in the tables of proportional parts had been given to 
tenths, as in Gauss's table, and the interpolated functions had been carried to 6 
places, besides the quantity xc' + (1 — x) c there would have been, arising from 
the neglect of the 7th place, a probable error of 0.025 ( m units of the 5th place) 
when the proportional part was given directly in the table, and of 0.035 if two 
entrances into the interpolation-table were made. This has so slight an effect 
that we pay no further attention to it, and assume that the error arising from 
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dropping the 6th place is 0.25. The last formula given above applies to this case, 
and for the average probable error we find 0.32. 

Class III. In Olney's or Davies's 6-place tables the difference for one second 
has a probable error of 0.0025 m units of the 6th place. This is multiplied by 
numbers ranging from o. 1 to 59.9. The average probable error of the product is 
thus 30 X 0.0025, or °-°7S- Combining this with the the probable error of the 
mantissa to which the product is added, and remembering that 0.25 is rejected 
when the result is limited to 6 places, we have 

Error = |/[(o.25) 2 + (0.075) 2 + (0.25) 2 ] = 0.36. 

The errors arising in the use of Schlomilch's 5-place table are nearly the 
same in magnitude as those of the tables in Class I. Let D be the apparent 
tabular difference and n the number of seconds. We really multiply Z>/6o (car- 
ried to hundredths) by n. But this is practically equivalent to multiplying «/6o 
(carried to hundredths) by D. Proceeding as under Case I, and including the 
error 0.075 which was lately mentioned, we find for the result 

Error = |/[(o.25) 2 + (o.8io) 2 (o.25) 2 + (0.075) 2 ] = 0.33. 

Hence functions interpolated from Schlomilch's table have a less probable error 
than those taken from Olney's or Davies's. This accords with the conclusion 
reached in Problem VI. 

PROBLEM XII. 

Given the logarithm of the tangent of an angle, to find the probable error of 
the corresponding interpolated angle. 

We divide the tables into the three classes given in Problem XI. 

Class I. The interpolation in this class of tables is like that in tables of 
logarithms of natural numbers. Hence the solution of this problem is similar to 
that of Problem II. If but one figure is interpolated we have 



*™->| [(#)' + («#]. 



If more than one figure is interpolated we have (see under Class I of Problem 
XI) instead of 1.98, 8.10/ 4, or 2.02. For this case 

Error = > ([f^. I c/] 2 + (o. 2 5) 2 ]. 

By substituting the values of D for different portions of the table, one can 
compute a series of values for the error. 

Class II. Let us assume that the proportional parts are given to tenths, as 
in Gauss's table. Suppose that the angles are carried accurately to thousandths 
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of a minute by computing d\D. Since the neglect of the ten-thousandths would 
introduce a probable error of only 60" X 0.00025, or 0V015 into the angle 
finally derived, and since the other errors are very much larger than this, we pay 
no further attention to it, but account the thousandths accurate. According to the 
reasoning in Problem II (remembering to use 0.810/4 instead of 0.793/4), we 
then have o.202/Z> as the probable error. But, in reducing the angle to seconds, 
this error would be multiplied by 60", giving the expression 12V12/Z?. The 
error thus far treated of arises from the inaccuracy of the tabular logarithms. 
But in using the tables of proportional parts given in Gauss, further error is 
introduced, because the proportional parts are carried to tenths only. Now, in 
general, if an angle be carried to hundredths of a second, three subtractions 
will be made, and at each subtraction a probable error of 0.025 °f a unit ' s intro- 
duced into the remainder. Hence the probable error of the final remainder is 
0.0251/3. Since the change of the proportional part for 1" is D/60, the effect 
in seconds of the error just mentioned is 0.025^/3 -i- Z>/6o, which equals 
1.51/3 ID = 2.6/D. The final value of the angle being restricted to seconds, a 
further error, whose probable value is o'.'25 is introduced. The three errors 
which we have discussed being independent of each other, the final result is 

Class III. In using Olney's or Davies's 6-place tables, the difference between 
the given logarithm and the tabular one has a probable error of 0.25. Denoting 
this difference by d, and the difference for 1" by D', since the latter has a prob- 
able error of 0.0025, we divide d ± 0.25 by D' ± 0.0025 and obtain the quotient 

d , 0.25 , . d 

D' + ~W ^ (°-°° 2 5) $n — ■ ■ ■ ■ 

We really adopt the value d\D' . Hence the error is 0.25 \D' =P 0.0025 d\D n . 
If the result is restricted to tenths of a second, as usual, another error whose 
probable value is 0.025 ' s introduced. The average value of djD 1 is 30, and the 
formula for the average probable error becomes 

or if D is taken as the difference for one minute, 

In employing Schlomilch's table, if the difference for 1" were carried to 4 or 5 
places, and the interpolated seconds to hundredths, the result would be the same 
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as if d were divided by D, the quotient carried to 4 or 5 places of decimals, and 
then multiplied by 60". 

But in this case the error due to c and c' would be, according to Problem II, 



I? n D ' 

whose probable value has been shown to be 0.202 \D when expressed as a 
decimal of a minute, or \2.\2\D, when expressed in seconds. As in Olney's 
table there is another error of 0.075 \D' due to tr, e inaccuracy of D' ; when the in- 
terpolated seconds are not carried beyond units, 0.25 is the probable value of 
the rejected figures. Thus the final formula is (1" being the unit) 

^[(^r+te ! r+ <->■]■ 

PROBLEM XIII. 

To find the maximum errors in Problem XI. 

Class I. The maximum error is 1.0 in the n' h place for «-place logarithms, as 
shown in Problem IV, 

Class II. For Houel's table no figure is rejected from the result given by 
interpolation. In addition to the error caused by the incorrectness of the tabular 
logarithms, one is introduced depending on the fact that each of the two func- 
tions taken from the table of proportional parts may have an error of 0.5 in the 
fifth place, and these errors may conspire. The formula then is 

Maximum error = xc' + (1 — x) c + 1.00 
= 1.50. 

If Gauss's table be used, the error due to cutting off the functions in the 
tables of proportional parts at tenths of a unit of the fifth decimal place, may 
amount to 0.05 -(- 0.05, or 0.10. Since all figures beyond the fifth place are re- 
jected, another error of 0.5 may arise. Hence we have 

Maximum error = xc' + (1 — x) c -f- 0.6 
= 1.10. 

Class III. In Olney's table an error arises from the neglect of second differ- 
ences, which is often quite appreciable, but will not be treated of. The difference 
for 1" may have an error of 0.005 m units of the sixth place. The product of 
this by 59.9 (the largest number of seconds ever employed) is 0.30. The maxi- 
mum error of the tabular mantissa to which the proportional part is applied is 
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O.50. The rejected figures beyond the sixth place may amount to 0.50. Hence 

Maximum error = 1.30. 

If Schlomilch's table be used, the error in the difference for 1", caused by 
terminating it at the second decimal place, may be 0.00^ of a unit of the fifth 
place. This multiplied by 59 gives 0.20. This table is also subject to the error 
of those of Class I. Therefore 

Maximum error = 1.20 

PROBLEM XIV. 

To find the maximum errors in Problem XII. 
Class I. From Problem V we get 

Maximum error = -i- . to 7 + 0.50, 

in which d is the minimum value of D, and / denotes the number of figures 
interpolated. The error is given in units of the last interpolated figure. For 
angles taken from the logarithmic tangents of Bremiker's 6-place table the for- 
mula gives 1.7 tenths of a second, or o'.'i7- 

Class II. We shall not discuss Hoiiel's table, because the value of the inter- 
polated seconds frequently depends upon the computer's judgment. 

In using Gauss's table let us suppose that the angles are carried to hun- 
dredths of a second, the process being to divide the difference d by the tabular 
difference D, and multiply the result by 60. Then by Problem II the error in sec- 
onds is 

6o"[^c' + l ^c-^*c(c> -c) +-^c'(c> -c)+ . . .]• 

Since the first two terms of the parenthesis are much larger than the others, we 
take them to find the condition for a maximum. Plainly, c' must equal c. Sub- 
stituting 0.50 for them, the foregoing expression becomes 



6o» [*£) , 



which is a maximum when D is a minimum. The least value of D in the tan- 
gent column is 25, and then the above expression = 1V20. But the employment 
of Gauss's tables of proportional parts causes another error, because the 
functions in them are carried to tenths only. A little experimentation shows 
that the maximum error due to this cause, when D is near 25, is o'.'20. In 
restricting the angle to the nearest second a possible error of o'.'5o is made, 
which added to those previously given makes I'.'tjo. It is hardly possible that 
this limit of inaccuracy is ever reached in practice. 
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Class III. In Olney's table we assume that the difference for i" has only the 
error 0.005, and that the next less tabular mantissa is 0.5 in error in the 6th 
place. Then Problem XII gives us 

Maximum error = — 'p- -f o"oo5 j, -s- 8' + o'.'os, 

in which 8' is the minimum value of £>', and the interpolated seconds are termi- 
nated at tenths. Since 8' = 4.21, the error is o"24. 

In employing Schlomilch's table there are three errors, the first due to the 
inaccuracy of the tabular logarithms, the second to the cutting off of D' at hun- 
dredths, and the third caused by restricting the interpolated angle to the nearest 
second. From previous discussions we readily get 

Maximum error = i"20 -f o'.'o03^ jA -+- 0' -f o'.'50. 

Assuming that d\8' = 57, and 8' = 0.42, this reduces to 2'.'i 5. In practice this 
error is never reached, and it is doubtful whether the limit 2'.'oo can be attained. 

problem xv. 

To compare the theories of Problems XI and XII with the facts. 

In testing the theories of these Problems, Bremiker's 6-place table was used 
as a specimen of those of Class I. 120 angles uniformly distributed between 5 
and 65 °, and carried to seconds, were chosen, and their logarithmic tangents 
taken from each of the five tables given below. The results were compared with 
those obtained from Schron's 7-place table, which was used as a standard. The 
averages of the errors, in units of the n' h place, for «-place logarithms are given 
below : — 

Actual Errors. Theoretical Errors. 

Bremiker, 0.27 0.31 

Houel, 0.39 0.39 

Gauss, 0.31 0.32 

Olney, 0.36 0.36 

Schlomilch, 0.31 0.33 

The average of the errors of the functions taken from Bremiker was un- 
usually small. 

In testing the theories of Problem XII, 120 logarithmic tangents were 
taken, and the corresponding angles, which ranged from 35 ° to 45°, were ob- 
tained. Schron was used again as the standard. For Bremiker the average 
value of D was 43, and for Gauss, 26. For Olney and for Schlomilch the aver- 
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age values of D' were 4.30 and 0.43 respectively. The average errors are tabu- 
lated below : — 

Actual Error. Theoretical Error. 

Bremiker, o'.'oso o'.'o53 

Gauss, 0.45 0.54 

Olney, 0.063 0.066 

Schlomilch, 0.50 0.56 

Being dissatisfied with the discrepancies in the cases of Gauss and Schlo- 
milch, I worked out another set of 120 for each, taking care to have the resulting 
seconds quite uniformly distributed throughout the minute. Gauss's error 
became 0V53, but Schlomilch's was 0V49 for forward interpolation and o'.'52 for 
backward interpolation. The seconds given by the interpolation backward dif- 
fered by a unit from those given by the forward, in 61 of the 120 angles; this 
portion of the table seems to be bewitched. 

Remark. — As it may be of interest to know the different times required in 
using various tables, I append those for obtaining the 1 20 angles from their loga- 
rithmic tangents. Schron's angles were carried to hundredths of a second, 
Bremiker's and Olney's to tenths, and those from the remaining tables to 
seconds. The degrees and minutes were written down in Schron's angles only, 
the seconds alone being recorded for the other tables. A multiplication table 
was used to facilitate interpolation from Olney. The seconds of each angle were 
found twice, once by forward interpolation, and once by backward, to insure 
their accuracy. 

First 40. Second 40. Third 40. Total. 

m. sec. m. sec. m. sec. h. m. sec. 
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Olney (6-place), 
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Gauss (5-place), 
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40 


Schlomilch (5-place), 
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10 





35 


20 


Gernerth (5-place), 
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i5 
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25 
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30 





14 


10 



In Gernerth's table the functions are given for each 10". 

Note. — Since logarithmic errors do not follow with exactness the laws of 
those treated of in the method of least squares, investigations based on the 
law assumed at the beginning of this series of articles cannot be expected to 
give results rigorously true; but their close approach to the truth, together with 
the difficulties of a rigid investigation of the subject, seem to justify the approx- 
imate method here employed. 



